Reciprocal regulatory balance within the CLEC16A-RNF41 mitophagy complex depends on an intrinsically disordered protein region.

CLEC16A is an E3 ubiquitin ligase that regulates mitochondrial quality control through mitophagy and is associated with over 20 human diseases. CLEC16A forms a complex with another E3 ligase, RNF41, and a ubiquitin-specific peptidase, USP8; however, regions that regulate CLEC16A activity or the assembly of the tripartite mitophagy regulatory complex are unknown. Here, we report that CLEC16A contains an internal intrinsically disordered protein region (IDPR) that is crucial for CLEC16A function and turnover. IDPRs lack a fixed secondary structure and possess emerging yet still equivocal roles in protein stability, interactions, and enzymatic activity. We find that the internal IDPR of CLEC16A is crucial for its degradation. CLEC16A turnover was promoted by RNF41, which binds and acts upon the internal IDPR to destabilize CLEC16A. Loss of this internal IDPR also destabilized the ubiquitin-dependent tripartite CLEC16A-RNF41-USP8 complex. Finally, the presence of an internal IDPR within CLEC16A was confirmed using NMR and CD spectroscopy. Together, our studies reveal that an IDPR is essential to control the reciprocal regulatory balance between CLEC16A and RNF41, which could be targeted to improve mitochondrial health in disease.

An intrinsically disordered protein region encoded by the human disease gene CLEC16A regulates mitophagy.

CLEC16A regulates mitochondrial health through mitophagy and is associated with over 20 human diseases. However, the key structural and functional regions of CLEC16A, and their relevance for human disease, remain unknown. Here, we report that a disease-associated CLEC16A variant lacks a C-terminal intrinsically disordered protein region (IDPR) that is critical for mitochondrial quality control. IDPRs comprise nearly half of the human proteome, yet their mechanistic roles in human disease are poorly understood. Using carbon detect NMR, we find that the CLEC16A C terminus lacks secondary structure, validating the presence of an IDPR. Loss of the CLEC16A C-terminal IDPR in vivo impairs mitophagy, mitochondrial function, and glucose-stimulated insulin secretion, ultimately causing glucose intolerance. Deletion of the CLEC16A C-terminal IDPR increases CLEC16A ubiquitination and degradation, thus impairing assembly of the mitophagy regulatory machinery. Importantly, CLEC16A stability is dependent on proline bias within the C-terminal IDPR, but not amino acid sequence order or charge. Together, we elucidate how an IDPR in CLEC16A regulates mitophagy and implicate pathogenic human gene variants that disrupt IDPRs as novel contributors to diabetes and other CLEC16A-associated diseases.Abbreviations : CAS: carbon-detect amino-acid specific; IDPR: intrinsically disordered protein region; MEFs: mouse embryonic fibroblasts; NMR: nuclear magnetic resonance.

Concepts and considerations for enhancing RNAi efficiency in phytopathogenic fungi for RNAi-based crop protection using nanocarrier-mediated dsRNA delivery systems.

Existing, emerging, and reemerging strains of phytopathogenic fungi pose a significant threat to agricultural productivity globally. This risk is further exacerbated by the lack of resistance source(s) in plants or a breakdown of resistance by pathogens through co-evolution. In recent years, attenuation of essential pathogen gene(s) via double-stranded (ds) RNA-mediated RNA interference (RNAi) in host plants, a phenomenon known as host-induced gene silencing, has gained significant attention as a way to combat pathogen attack. Yet, due to biosafety concerns regarding transgenics, country-specific GMO legislation has limited the practical application of desirable attributes in plants. The topical application of dsRNA/siRNA targeting essential fungal gene(s) through spray-induced gene silencing (SIGS) on host plants has opened up a transgene-free avenue for crop protection. However, several factors influence the outcome of RNAi, including but not limited to RNAi mechanism in plant/fungi, dsRNA/siRNA uptake efficiency, dsRNA/siRNA design parameters, dsRNA stability and delivery strategy, off-target effects, etc. This review emphasizes the significance of these factors and suggests appropriate measures to consider while designing in silico and in vitro experiments for successful RNAi in open-field conditions. We also highlight prospective nanoparticles as smart delivery vehicles for deploying RNAi molecules in plant systems for long-term crop protection and ecosystem compatibility. Lastly, we provide specific directions for future investigations that focus on blending nanotechnology and RNAi-based fungal control for practical applications.

Distinct conformational dynamics and allosteric networks in alpha tryptophan synthase during active catalysis.

Experimental observations of enzymes under active turnover conditions have brought new insight into the role of protein motions and allosteric networks in catalysis. Many of these studies characterize enzymes under dynamic chemical equilibrium conditions, in which the enzyme is actively catalyzing both the forward and reverse reactions during data acquisition. We have previously analyzed conformational dynamics and allosteric networks of the alpha subunit of tryptophan synthase under such conditions using NMR. We have proposed that this working state represents a four to one ratio of the enzyme bound with the indole-3-glycerol phosphate substrate (E:IGP) to the enzyme bound with the products indole and glyceraldehyde-3-phosphate (E:indole:G3P). Here, we analyze the inactive D60N variant to deconvolute the contributions of the substrate- and products-bound states to the working state. While the D60N substitution itself induces small structural and dynamic changes, the D60N E:IGP and E:indole:G3P states cannot entirely account for the conformational dynamics and allosteric networks present in the working state. The act of chemical bond breakage and/or formation, or possibly the generation of an intermediate, may alter the structure and dynamics present in the working state. As the enzyme transitions from the substrate-bound to the products-bound state, millisecond conformational exchange processes are quenched and new allosteric connections are made between the alpha active site and the surface which interfaces with the beta subunit. The structural ordering of the enzyme and these new allosteric connections may be important in coordinating the channeling of the indole product into the beta subunit.

Different Solvent and Conformational Entropy Contributions to the Allosteric Activation and Inhibition Mechanisms of Yeast Chorismate Mutase.

Allosteric regulation is important in many biological processes, including cell signaling, gene regulation, and metabolism. Saccharomyces cerevisiae chorismate mutase (ScCM) is a key homodimeric enzyme in the shikimate pathway responsible for the generation of aromatic amino acids, where it is allosterically inhibited and activated by Tyr and Trp, respectively. Our previous studies indicated that binding of both allosteric effectors is negatively cooperative, that is binding at one allosteric binding site discourages binding at the other, due to the entropic penalty of binding the second allosteric effector. We utilized variable temperature isothermal titration calorimetry (ITC) and nuclear magnetic resonance (NMR) experiments to better understand the entropic contributions to allosteric effector binding, including changes to solvent entropy and protein conformational entropy. Upon binding either Tyr or Trp, ScCM experiences a quenching of motions on the picosecond-to-nanosecond time scale, which we could relate to a loss of protein conformational entropy. Further ITC and NMR studies were consistent with the Tyr-bound form of ScCM being associated with more water molecules compared to the Trp-bound form and Tyr binding being associated with a less positive solvent entropy change. These studies provide insight into the role of structural dynamics in ScCM function and highlight the importance of solvent entropy changes in allosteric regulation, a historically underappreciated concept.

Diastereoselective multi-component tandem condensation: synthesis of 2-amino-4-(2-furanone)-4H-chromene-3-carbonitriles.

A general strategy for a one-pot stereoselective synthesis of 2-amino-4-(2-furanone)-4H-chromene-3-carbonitriles by reaction of salicylaldehyde, malononitrile and butenolides via a tandem Knoevenagel/Pinner/vinylogous Michael condensation is presented. The ß,γ-butenolides gave a syn-selective MCR adduct with a dr up to 11.5 : 1. The mechanistic insight into the MCR was obtained by DFT calculations.

Coordinated Network Changes across the Catalytic Cycle of Alpha Tryptophan Synthase.

Networks of noncovalent interactions are important for protein structural dynamics. We used nuclear magnetic resonance chemical shift covariance analyses on an inactive variant of the alpha subunit of tryptophan synthase to map amino acid interaction networks across its catalytic cycle. Although some network connections were common to every enzyme state, many of the network connections strengthened or weakened over the catalytic cycle; these changes were highly coordinated. These results suggest a higher level of network organization. Our analyses identified periodic, second-order networks that show highly coordinated interaction changes across the catalytic cycle. These periodic networks may help synchronize the sequence of structural transitions necessary for enzyme function. Molecular dynamics simulations identified interaction changes across the catalytic cycle, including those involving the catalytic residue Glu49, which may help drive other interaction changes throughout the enzyme structure. Similar periodic networks may direct structural transitions and allosteric interactions in other proteins.

Solution Ensemble of the C-Terminal Domain from the Transcription Factor Pdx1 Resembles an Excluded Volume Polymer.

The pancreatic and duodenal homeobox 1 (Pdx1) is an essential pancreatic transcription factor. The C-terminal intrinsically disordered domain of Pdx1 (Pdx1-C) has a heavily biased amino acid composition; most notably, 18 of 83 residues are proline, including a hexaproline cluster near the middle of the chain. For these reasons, Pdx1-C is an attractive target for structure characterization, given the availability of suitable methods. To determine the solution ensembles of disordered proteins, we have developed a suite of 13C direct-detect NMR experiments that provide high spectral quality, even in the presence of strong proline enrichment. Here, we have extended our suite of NMR experiments to include four new pulse programs designed to record backbone residual dipolar couplings in a 13C,15N-CON detection format. Using our NMR strategy, in combination with small-angle X-ray scattering measurements and Monte Carlo simulations, we have determined that Pdx1-C is extended in solution, with a radius of gyration and internal scaling similar to that of an excluded volume polymer, and a subtle tendency toward a collapsed structure to the N-terminal side of the hexaproline sequence. This structure leaves Pdx1-C exposed for interactions with trans-regulatory co-factors that contribute with Pdx1 to transcription control in the cell.

Millisecond Timescale Motions Connect Amino Acid Interaction Networks in Alpha Tryptophan Synthase.

Tryptophan synthase is a model system for understanding allosteric regulation within enzyme complexes. Amino acid interaction networks were previously delineated in the isolated alpha subunit (αTS) in the absence of the beta subunit (ßTS). The amino acid interaction networks were different between the ligand-free enzyme and the enzyme actively catalyzing turnover. Previous X-ray crystallography studies indicated only minor localized changes when ligands bind αTS, and so, structural changes alone could not explain the changes to the amino acid interaction networks. We hypothesized that the network changes could instead be related to changes in conformational dynamics. As such, we conducted nuclear magnetic resonance relaxation studies on different substrate- and products-bound complexes of αTS. Specifically, we collected 15N R2 relaxation dispersion data that reports on microsecond-to-millisecond timescale motion of backbone amide groups. These experiments indicated that there are conformational exchange events throughout αTS. Substrate and product binding change specific motional pathways throughout the enzyme, and these pathways connect the previously identified network residues. These pathways reach the αTS/ßTS binding interface, suggesting that the identified dynamic networks may also be important for communication with the ßTS subunit.

Assigning methyl resonances for protein solution-state NMR studies.

Solution-state NMR is an important tool for studying protein structure and function. The ability to probe methyl groups has substantially expanded the scope of proteins accessible by NMR spectroscopy, including facilitating study of proteins and complexes greater than 100 kDa in size. While the toolset for studying protein structure and dynamics by NMR continues to grow, a major rate-limiting step in these studies is the initial resonance assignments, especially for larger (>50 kDa) proteins. In this practical review, we present strategies to efficiently isotopically label proteins, delineate NMR pulse sequences that can be used to determine methyl resonance assignments in the presence and absence of backbone assignments, and outline computational methods for NMR data analysis. We use our experiences from assigning methyl resonances for the aromatic biosynthetic enzymes tryptophan synthase and chorismate mutase to provide advice for all stages of experimental set-up and data analysis.

Discrete-State Kinetics Model for NMR-Based Analysis of Protein Translocation on DNA at Equilibrium.

In the target DNA search process, sequence-specific DNA-binding proteins first nonspecifically bind to DNA and stochastically move from one site to another before reaching their targets. To rigorously assess how the translocation process influences NMR signals from proteins interacting with nonspecific DNA, we incorporated a discrete-state kinetic model for protein translocation on DNA into the McConnell equation. Using this equation, we simulated line shapes of NMR signals from proteins undergoing translocations on DNA through sliding, dissociation/reassociation, and intersegment transfer. Through this analysis, we validated an existing NMR approach for kinetic investigations of protein translocation on DNA, which utilizes NMR line shapes of two nonspecific DNA-protein complexes and their mixture. We found that, despite its use of simplistic two-state approximation neglecting the presence of many microscopic states, the previously proposed NMR approach provides accurate kinetic information on the intermolecular translocations of proteins between two DNA molecules. Interestingly, our results suggest that the same NMR approach can also provide qualitative information about the one-dimensional diffusion coefficient for proteins sliding on DNA.

Engineering double-stranded RNA binding activity into the Drosha double-stranded RNA binding domain results in a loss of microRNA processing function.

Canonical processing of miRNA begins in the nucleus with the Microprocessor complex, which is minimally composed of the RNase III enzyme Drosha and two copies of its cofactor protein DGCR8. In structural analogy to most RNase III enzymes, Drosha possesses a modular domain with the double-stranded RNA binding domain (dsRBD) fold. Unlike the dsRBDs found in most members of the RNase III family, the Drosha-dsRBD does not display double-stranded RNA binding activity; perhaps related to this, the Drosha-dsRBD amino acid sequence does not conform well to the canonical patterns expected for a dsRBD. In this article, we investigate the impact on miRNA processing of engineering double-stranded RNA binding activity into Drosha's non-canonical dsRBD. Our findings corroborate previous studies that have demonstrated the Drosha-dsRBD is necessary for miRNA processing and suggest that the amino acid composition in the second α-helix of the domain is critical to support its evolved function.

Assessing Coupled Protein Folding and Binding Through Temperature-Dependent Isothermal Titration Calorimetry.

Broad interest in the thermodynamic driving forces of coupled macromolecular folding and binding is motivated by the prevalence of disorder-to-order transitions observed when intrinsically disordered proteins (IDPs) bind to their partners. Isothermal titration calorimetry (ITC) is one of the few methods available for completely evaluating the thermodynamic parameters describing a protein-ligand binding event. Significantly, when the effective ΔH° for the coupled folding and binding process is determined by ITC in a temperature series, the constant-pressure heat capacity change (ΔCp) associated with these coupled equilibria is experimentally accessible, offering a unique opportunity to investigate the driving forces behind them. Notably, each of these molecular-scale events is often accompanied by strongly temperature-dependent enthalpy changes, even over the narrow temperature range experimentally accessible for biomolecules, making single temperature determinations of ΔH° less informative than typically assumed. Here, we will document the procedures we have adopted in our laboratory for designing, executing, and globally analyzing temperature-dependent ITC studies of coupled folding and binding in IDP interactions. As a biologically significant example, our recent evaluation of temperature-dependent interactions between the disordered tail of FCP1 and the winged-helix domain from Rap74 will be presented. Emphasis will be placed on the use of publically available analysis programs written in MATLAB that facilitate quantification of the thermodynamic forces governing IDP interactions. Although motivated from the perspective of IDPs, the experimental design principles and data fitting procedures presented here are general to the study of most noncooperative ligand binding equilibria.

Generating NMR chemical shift assignments of intrinsically disordered proteins using carbon-detected NMR methods.

There is an extraordinary need to describe the structures of intrinsically disordered proteins (IDPs) due to their role in various biological processes involved in signaling and transcription. However, general study of IDPs by NMR spectroscopy is limited by the poor (1)H amide chemical shift dispersion typically observed in their spectra. Recently, (13)C direct-detected NMR spectroscopy has been recognized as enabling broad structural study of IDPs. Most notably, multidimensional experiments based on the (15)N,(13)C CON spectrum make complete chemical shift assignment feasible. Here we document a collection of NMR-based tools that efficiently lead to chemical shift assignment of IDPs, motivated by a case study of the C-terminal disordered region from the human pancreatic transcription factor Pdx1. Our strategy builds on the combination of two three-dimensional (3D) experiments, (HN-flip)N(CA)CON and 3D (HN-flip)N(CA)NCO, that enable daisy chain connections to be built along the IDP backbone, facilitated by acquisition of amino acid-specific (15)N,(13)C CON-detected experiments. Assignments are completed through carbon-detected, total correlation spectroscopy (TOCSY)-based side chain chemical shift measurement. Conducting our study required producing valuable modifications to many previously published pulse sequences, motivating us to announce the creation of a database of our pulse programs, which we make freely available through our website.

Real-time kinetics of high-mobility group box 1 (HMGB1) oxidation in extracellular fluids studied by in situ protein NMR spectroscopy.

Some extracellular proteins are initially secreted in reduced forms via a non-canonical pathway bypassing the endoplasmic reticulum and become oxidized in the extracellular space. One such protein is HMGB1 (high-mobility group box 1). Extracellular HMGB1 has different redox states that play distinct roles in inflammation. Using a unique NMR-based approach, we have investigated the kinetics of HMGB1 oxidation and the half-lives of all-thiol and disulfide HMGB1 species in serum, saliva, and cell culture medium. In this approach, salt-free lyophilized (15)N-labeled all-thiol HMGB1 was dissolved in actual extracellular fluids, and the oxidation and clearance kinetics were monitored in situ by recording a series of heteronuclear (1)H-(15)N correlation spectra. We found that the half-life depends significantly on the extracellular environment. For example, the half-life of all-thiol HMGB1 ranged from ~17 min (in human serum and saliva) to 3 h (in prostate cancer cell culture medium). Furthermore, the binding of ligands (glycyrrhizin and heparin) to HMGB1 significantly modulated the oxidation kinetics. Thus, the balance between the roles of all-thiol and disulfide HMGB1 proteins depends significantly on the extracellular environment and can also be artificially modulated by ligands. This is important because extracellular HMGB1 has been suggested as a therapeutic target for inflammatory diseases and cancer. Our work demonstrates that the in situ protein NMR approach is powerful for investigating the behavior of proteins in actual extracellular fluids containing an enormous number of different molecules.

Ensemble analysis of primary microRNA structure reveals an extensive capacity to deform near the Drosha cleavage site.

Most noncoding RNAs function properly only when folded into complex three-dimensional (3D) structures, but the experimental determination of these structures remains challenging. Understanding of primary microRNA (miRNA) maturation is currently limited by a lack of determined structures for nonprocessed forms of the RNA. SHAPE chemistry efficiently determines RNA secondary structural information with single-nucleotide resolution, providing constraints suitable for input into MC-Pipeline for refinement of 3D structure models. Here we combine these approaches to analyze three structurally diverse primary microRNAs, revealing deviations from canonical double-stranded RNA structure in the stem adjacent to the Drosha cut site for all three. The necessity of these deformable sites for efficient processing is demonstrated through Drosha processing assays. The structure models generated herein support the hypothesis that deformable sequences spaced roughly once per turn of A-form helix, created by noncanonical structure elements, combine with the necessary single-stranded RNA-double-stranded RNA junction to define the correct Drosha cleavage site.

The role of human Dicer-dsRBD in processing small regulatory RNAs.

One of the most exciting recent developments in RNA biology has been the discovery of small non-coding RNAs that affect gene expression through the RNA interference (RNAi) mechanism. Two major classes of RNAs involved in RNAi are small interfering RNA (siRNA) and microRNA (miRNA). Dicer, an RNase III enzyme, plays a central role in the RNAi pathway by cleaving precursors of both of these classes of RNAs to form mature siRNAs and miRNAs, which are then loaded into the RNA-induced silencing complex (RISC). miRNA and siRNA precursors are quite structurally distinct; miRNA precursors are short, imperfect hairpins while siRNA precursors are long, perfect duplexes. Nonetheless, Dicer is able to process both. Dicer, like the majority of RNase III enzymes, contains a dsRNA binding domain (dsRBD), but the data are sparse on the exact role this domain plays in the mechanism of Dicer binding and cleavage. To further explore the role of human Dicer-dsRBD in the RNAi pathway, we determined its binding affinity to various RNAs modeling both miRNA and siRNA precursors. Our study shows that Dicer-dsRBD is an avid binder of dsRNA, but its binding is only minimally influenced by a single-stranded - double-stranded junction caused by large terminal loops observed in miRNA precursors. Thus, the Dicer-dsRBD contributes directly to substrate binding but not to the mechanism of differentiating between pre-miRNA and pre-siRNA. In addition, NMR spin relaxation and MD simulations provide an overview of the role that dynamics contribute to the binding mechanism. We compare this current study with our previous studies of the dsRBDs from Drosha and DGCR8 to give a dynamic profile of dsRBDs in their apo-state and a mechanistic view of dsRNA binding by dsRBDs in general.

Asymmetrical roles of zinc fingers in dynamic DNA-scanning process by the inducible transcription factor Egr-1.

Egr-1 is an inducible transcription factor that recognizes 9-bp target DNA sites via three zinc finger domains and activates genes in response to cellular stimuli such as synaptic signals and vascular stresses. Using spectroscopic and computational approaches, we have studied structural, dynamic, and kinetic aspects of the DNA-scanning process in which Egr-1 is nonspecifically bound to DNA and perpetually changes its location on DNA. Our NMR data indicate that Egr-1 undergoes highly dynamic domain motions when scanning DNA. In particular, the zinc finger 1 (ZF1) of Egr-1 in the nonspecific complex is mainly dissociated from DNA and undergoes collective motions on a nanosecond timescale, whereas zinc fingers 2 and 3 (ZF2 and ZF3, respectively) are bound to DNA. This was totally unexpected because the previous crystallographic studies of the specific complex indicated that all of Egr-1's three zinc fingers are equally involved in binding to a target DNA site. Mutations that are expected to enhance ZF1's interactions with DNA and with ZF2 were found to reduce ZF1's domain motions in the nonspecific complex suggesting that these interactions dictate the dynamic behavior of ZF1. By experiment and computation, we have also investigated kinetics of Egr-1's translocation between two nonspecific DNA duplexes. Our data on the wild type and mutant proteins suggest that the domain dynamics facilitate Egr-1's intersegment transfer that involves transient bridging of two DNA sites. These results shed light on asymmetrical roles of the zinc finger domains for Egr-1 to scan DNA efficiently in the nucleus.

NMR studies of translocation of the Zif268 protein between its target DNA Sites.

Zif268 is a zinc-finger protein containing three Cys(2)-His(2)-type zinc-finger domains that bind the target DNA sequence GCGTGGGCG in a cooperative manner. In this work, we characterized translocation of the Zif268 protein between its target DNA sites using NMR spectroscopy. The residual dipolar coupling data and NMR chemical shift data suggested that the structure of the sequence-specific complex between Zif268 and its target DNA in solution is the same as the crystal structure. Using two-dimensional heteronuclear (1)H-(15)N correlation spectra recorded with the fast acquisition method, we analyzed the kinetics of the process in which the Zif268 protein transfers from a target site to another on a different DNA molecule on a minute to hour time scale. By globally fitting the time-course data collected at some different DNA concentrations, we determined the dissociation rate constant for the specific complex and the second-order rate constant for direct transfer of Zif268 from one target site to another. Interestingly, direct transfer of the Zif268 protein between its target sites is >30000-fold slower than corresponding direct transfers of the HoxD9 and the Oct-1 proteins, although the affinities of the three proteins to their target DNA sites are comparable. We also analyzed translocation of the Zif268 protein between two target sites on the same DNA molecules. The populations of the proteins bound to the target sites were found to depend on locations and orientations of the target sites.

Redox properties of the A-domain of the HMGB1 protein.

The High Mobility Group B1 (HMGB1) protein plays important roles in both intracellular (reductive) and extracellular (oxidative) environments. We have carried out quantitative investigations of the redox chemistry involving Cys22 and Cys44 of the HMGB1 A-domain, which form an intramolecular disulfide bond. Using NMR spectroscopy, we analyzed the real-time kinetics of the redox reactions for the A-domain in glutathione and thioredoxin systems, and also determined the standard redox potential. Thermodynamic experiments showed that the Cys22-Cys44 disulfide bond stabilizes the folded state of the protein. These data suggest that the oxidized HMGB1 may accumulate even in cells under oxidative stress.